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SUMMARY 

Hughes Research Laboratories has fabricated state-of-the-art planar 
star couplers for JPL for distributing optical radiation between optical 
fibers. These devices were fabricated by means of a new ion-exchange process. 
An 8x8 planar transmission star coupler has been packaged for evaluation pur- 
poses, with sixteen fiber connectors and sixteen pigtails. An additional 
transmission star coupler and an eight-port reflection star coupler have 
been fabricated and delivered; eight-fiber ribbons are rigidly attached to 
these couplers. A planar coupler with silicon guides and a parallel channel 
guide with pigtails have been fabricated and delivered. Optical measurements 
of the transmissior star couplers have been made and are 
included in this final report. 


vii 


INTRODUCTION 

Low optical attenuation has made the optical fiber waveguide the leading 
contender as a transmission medium for many future communication and data 
handling systems. Low loss fibers have significant advantages over metallic 
conductors. Including large bandwidth, small size and weight. Immunity from 
crosstalk, ground loops and electromagnetic Interference. Optical fibers 
thus appear to be an ideal medium for communication and data links operating 
in harsh environments such as experienced by spacecraft. 

Development of specific components are required to take full advantage 
of optical fiber links. Such a component Is the optical data distribution 
device, or coupler, which partitions radiation contained in one fiber to 
other fibers. This contract was directed towards fabricating state-of-the-art 
couplers for performing the data distribution task. We have fabricated a 
number of such couplers and have delivered the best planar couplers that we 
have fabricated to date on this contract. The design of the coupler, the 
chemical process used, the assembly process, and the optical measurement 
made on the delivered items are described here. 
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PLANAR STAR COUPLER DESIGN 

The design of the transmission star coupler delivered under this contract 
Is shown in Figure 1. The planar star Is shown together with the Inter- 
connecting ribbons of eight fibers. The planar star consists of the inter- 
mediate channel waveguides, expansion and collecting horns, and the mixing 
region. These areas are processed by Ion exchange and have a higher refrac- 
tive index than the adjoining substrate, thereby permitting light guidance in 
these regions. 

The geometric design for equal distribution of light from one fiber into 
the fibers at the other end of the star was determined by a computer ray 
tracing program. The computer program is described in detail in Appendix A. 

A step refractive index profile was assumed for the waveguide structure. As 
many as 10,000 rays were launched and tracked through the planar structure 
during the computer program. The computer program determined the number of 
rtys partitioned into each exit port as the parameters (horn angle, inter- 
mediate channel width and spacing, mixing length) of the planar structure were 
varied. The program also summed the number of rays lost from the waveguide. 

The calculations and practical considerations of fiber size resulted in 
a planar transmission star with the following dimensions: 

Length of star, 10cm, 

Width of star, 0.2 cm, 

Intermediate channel width (mask opening), 47 urn. 

Separation of channel, center-to-center, 24? pm. 

Angle of horn, 2°, 

Length of mixing region, 8 cm. 
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Fig, 1, Eight-port transmission star coupler design 


The star structure was Ion-exchanged Into a plate glass substrate. 
Although the step refractive Index assumption Is not realized In the Ion 
exchange process, the ray tracing computer calculation (with step Index 
profile) served as a guide In the design of the star. With the dimensions 
given above, the computer result with step Index profile predicts a nearly 
equal partitioning of rays Into the exit ports with a loss of ^ 20 %. The 
observed partitioning In a fabricated device, however, deviates from this. 

It then became clear that computer programs must consider the actual index 
profile. To this end, we initiated a ray tracing calculation (under Hughes 
IR&D) for a parallel channel with a graded index profile, where the input 
and output fibers also have index profile. The results are encouraging. 

The calculated throughput of the channel appears to correlate closely with 
the measured throughput. 

A brief description of the ray tracing program for the parallel channel 
graded index structure is described here. We have computed the throughput 
power due to coupling a graded Index fiber to graded parallel channel to 
graded index fiber. This theoretical result serves as a guideline for evalu- 
ating planar couplers fabricate^ by the single ion exchange process. 

The assumptions and procedures used in the computer program are: 

The ion exchanged waveguide has been modeled as a stepwise graded 
index structure consisting of 20 approximately el 1 iptical ly shaped 
layers, each layer differing in refractive index. The refractive 
index varies linearly with layer depth from 1.532 at the ion exchange 
interface to 1.523 at the deepest layer. The length of the minor axis 
of an elliptical layer is taken to be twice the channel depth of that 
layer. The length of the major axis of each of the elliptical layer 
is set equal to the opening of the photomask (40 urn) plus the length 
of the mirror axis. The largest ellipse has a minor axis equal to 
twice the maximum waveguide depth. 

The graded index fiber to which we couple the planar structure is assumed 
to have the following properties, as a source of rays to the planar 
structure: 
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1) Ray density over the core surface Is uniform, 

2) The maximum angular spread of the rays varies from the maximum 
guide numerical aperture at the center of the core NA(0) to zero 
spread at the core edge. In accordance with the formula 

NA{r) - (1 - ar 2 ) NA(0) 

where r 1$ the radial distance from the center of the core. The 
fiber NA(0) Is assumed to be 0.2. 

The results of the calculations are shown In Figure 2 where the coupled 
output power is graphed as a function of the relative position of the fiber 
to the planar channel. The upper two curves describe the fractional power at 
the exit (right) end of the channel due to radiation coupled from the fiber. 

The dotted line represents the output for a channel depth of three core radii, 
and the continuous line the output for a depth of two core radii. The maxima 
occurs when the ratio of the core center to surface distance f core radius is 
unity, l.e. , when the upper surface of the fiber core is at the same level as 
the channel surface. 

The lower two curves of Figure 2 show the power coupled from the graded 
Index fiber to the channel and then to the output fiber to the right. The 
upper continuous curve shows the output when the planar channel is two core 
radii deep, the dotted curve the output for a channel three core radii deep. 

The results suggest that the channel depth should be less than or equal to 
the core diameter for maximum power throughput. The maxima of the upper curve 
occurs when the upper surface of the core Is slightly above the planar surface; 
approximately 5555 of the Input light should then be captured in the modes of 
the exit fiber. 

It is clear that ray tracing through a structure such as the transmission 
star requires a somewhat more complex program than that for the parallel channel. 
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RATIO: CORE CENTER TO SURFACE DISTANCE * CORE RADIUS 


Fig. 2. Calculated power throughput for (a) flber-to-channel 
coupling and (b) flber-to-channel -to-flber coupling* 
using graded Index parallel channel and graded Index 
fiber. Fiber radius * 31.25 um; channel An * 0.009; 
graded index fiber NA * 0.2; 40 um mask openlnci. 
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Nevertheless, we expect to address this problem In the Immediate future In 
order to be able to correlate with the experimental results and to design 
Improved stars and other planar devices. 
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PLANAR COUPLER PROCESSING 

The Ion-exchange process for fabricating the planar devices are described 
here. Elqht-port transmission star couplers were fabricated by the process 
described In Appendix B. Aluminum films were evaporated on microscope or plate 
glass slides, then photoetched with open area corresponding to the planar 
star design described above. The slide was then dipped in a molten solution 
of L^SO^-I^SO^ at 585°C for twenty (20) minutes to ion-exchange Na of the 
glass with Li of the melt. The slides were then sawed, ground and polished. 

We have fabricated over fifty experimental star couplers using this 
procedure. Most of these were fabricated to test the ion-exchange exposure 
(temperature and time) and to test the photolithographic process. The earlier 
stars exhibited considerable scattering and loss. We examined the star coupler 
fabrication process in further detail to determine the sources of loss. The 
result has been an improvement in the quality of the star couplers. First, 
we determined that, for earlier fabrications, the temperature controller for 
the Li^SO^-I^SO^ melt was not properly calibrated, yielding a 20°C higher 
temperature than Indicated. The higher temperature tends to corrode the A1 
masks more readily resulting In poorly formed guides. The temperature con- 
troller was recalibrated. Secondly, we have changed from negative to positive 
photoresists, using a dark field mask. The photoresist Is now dipcoated, 
resulting in more uniform layer of photoresist on the aluminum and therefore 
more uniform etching. (Previously, spin coating of resist on the long samples 
resulted in a wavy, non-uniform photoresist layer.) The result has been a 
clean aluminum mask formation on the qlass. More Importantly, the use of 
positive photoresist process has eliminated the polymerized photoresist which 
remains on the aluminum surface with the negative resist process. It is 
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possible that during previous processing the polymerized photoresist had 
decomposed Into products obscuring the star channels and resulting In incomplete 
ion exchange. 

We fabricated two sets of ten 8x8 transmission star couplers, one set 
with negative resist and one set with positive resist. Various temperatures 
and ion exchange immersion times were used. The results were dramatically 
different. The negative resist couplers exhibited poor light transmission. 

The positive resist process haS resulted in the best star coupler fabricated 
to date with throughput varying from better than -9 dB to -11 dB. 


FIBER RIBBON ASSEMBLY 

Two ribbons of eighteen (18) fibers (including fibers as spacers) to 
mate with the 8-channel star were fabricated. The cross section of one of 
the ribbons Is shown In Figure 3a. Assuming that the fiber diameters were 
constant at 125 pm, an eight-port star mask was specified with interchannel 
spacing of 250 pm. Figure 3b Illustrates dimensional mismatch arising from 
variations in fiber diameter and improper mask dimensions; the fiber spacings 
are greater than the channel spacings. Due to fluctuations in fiber diameter, 
the above ribbon fabrication procedure has been discontinued and replaced 
with a new process. We now use etched silicon alignment grooves to fabricate 
the fiber ribbon. 

We successfully fabricated several ribbons of eight fibers for mating 
with the 8x8 planar star couplers. Etched silicon channels or spacers were 
fabricated by a process used in semiconductor technology. Etching of Si in 
the 100 direction proceeds at a fast rate relative to other directions. 

Channel masks with openings to match the dimensions of the star mask channels 
were fabricated. Using these masks, deep V grooves were etched into Si for 
aligning the fibers. Eight fibers were then Inserted Into the etched silicon 
channels and epoxied into place; the end of the unit was then ground and 
polished. Figure 4 shows the end view of the ribbon of eight fibers. 

The ribbons of eight fibers were then carefully positioned at the ends 
of the eight-port transmission star coupler for maximum light throughput. 

The ribbons were then epoxied into place to form a rigid unit, as shown in 
Fig. 5a, b. One such unit was then boxed with connectors at the Hughes 
Connector Division in Irvine, Figures 6 and 7. (This unit is a prototype 
[or brassboard] and does not suggest any final packaging design.) 
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Fin. 3a. End view of fiber ribbon connector. 
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Fig. 3b. Photograph of the input end of a star coupler 

with mating fiber ribbon connector. Mote align- 
ment mismatch. 
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Fig. 4 . 


End view of ribbon of light fibers epoxied in etched 
silicon spacers. 
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Fig. 5a. 


Planar transmission star coupler with eight-fiber 
ribbon epoxied to each end. 
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Fig. 5b. Closeup of transmission star fiber ribbon connection. 
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Fig. 6. Planar transmission star in box with fiber connectors 
(cover removed). 
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OPTICAL MEASUREMENTS 

Optical measurements of the refractive Index profile of the Ion- 
exchanged channel, flber-to-channel coupling effects, and the transmission 
star coupler throughputs are described here. 

Refractive Index Profile 

The refractive Index profile of the Ion-exchanged channel was examined. 

An Interference microscope photograph of a thin section of three channels 
is shown In Figure 8. The highest Index Is 1.532, the region that directly 
contacted the melt, and the lowest Is 1.523, that of the glass substrate. 

For a step Index profile, these indices would correspond to an NA * .16. 

Firu^e 8 clearly illustrates that Ion exchange occurs under the A1 -masked 
area, diffusing laterally as well as normally Into the glass. The diffusion 
distance normal to the surface is ^100 pm. Smaller mask openings would result 
in a more circular refractive index pattern. 

We conclude that An = 0.009 is sufficiently large to be able to match 
the An of a common commercial fiber. The spatial or cross-sectional mismatch 
of indices between the fiber and the channel will, however, decrease the 
coupling efficiency to some extent. Processes for attaining closer cross- 
sectional index matching are under investigation at HRL. 

Optimization of the Flber-to-Channel Coupling 

In order to determine the optimum width of a constant width channel for 
best fiber-to-channel-to-flber and to determine the flber-to-channel coupling 
losses, a series of measurements were made on existing slides of straight 
channels of varying widths. Channel widths (of the mask) of 25, 40, 50, 

65, 80, 90, 250, 325, 500, and 620 pm and channel lengths of one inch were 
available on previously fabricated slides. The slides had been dipped into 


19 





Fig. 8. Interference microscope photograph of a star 
coupler cross section showing change of index 
of refraction. Center-to-center channel spacina 
and mask opening width are shown alonq with the 
indices of the substrate and ion exchanged channels. 
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the melt for twenty (20) minutes. The Ion exchanged regions for various 
channel widths of the mask are shown schematically In the Inset of Fig. 9. 

The diffusion of LI ions underneath the masks is Indicated. The diffusion 
depth (MOO pm) Is determined by length of dipping (20 minutes) and melt 
temperature (580°C). Graded Index fibers of 63 pm core and 125 pm o.d. were 
butt coupled to each end of a channel. 

Our measurements Indicated that maximum throughput for fiber to channel 
to fiber (M5%) occurred over a broad range of widths, 45-80 urn, shown in 
Fig. 9 and indicated by F/C/F. The flber-to-channel output versus channel 
widths are indicated by F/C In Fig. 9; the near field is focussed on a detector 
with a lens, with an iris to isolate the channel throughput from scattered 
light. The efficiency of coupling from the output channel to the fiber can 
then be estimated by taking the ratio of (F/C/F) to (F/C); the ratio is 
shown in Fig. 10 as a function of channel width, with a maximum of \75 f/ . 

In summary, we conclude that: 

t Fiber-to-channel-to-flber throughputs were measured to have a broad 
maxima of M5% for channel widths of 45 to 80 nm for the fabricated 
channels. This result includes channel scattering losses and channel - 
to-fiber coupling losses resulting from refractive index profile 
mismatch at both the Input and output. The graded index fiber has 
core diameter of 63 pm, o.d. of 130 um, and an NA of 0.14. Closer 
matching of profiles are desired, such as with buried channels. 

t The output channel -to-fiber losses for channel widths of 45 to 80 nm 
are »-25 

• T rom the straight channel experiments, we can estimate that the 
power entering a star coupler is >577. with present processing. We 
expect that a well fabricated star coupler will have total throughput 
somewhat greater than the above value. 
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Fig. 9. Power transmission ratios of constant width channel vs. 
channel widths for flber-to-channel coupling (F/C) and 
flber-to-channe 1 to-flber coupling (F/C/F). 
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Fig. 10. Channel -to- fiber coupling efficiency versus channel widths. 
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Transmission Star Throughput 

As mentioned above, the last set of ten transmission star couplers 
fabricated with the positive resist process has yielded the best quality 
stars thus far, with relatively low surface scattering. Plates containing 
the stars were Immersed In the LlgSO^-KgSO^ melt for the following temperature 
and time Intervals: 605®C, 5 min; 605°C, 10 min.; 605°C, 15 min.; 585°C, 

20 min.; 585°C, 40 min. Two stars were located on each plate. The plate 
Immersed for 20 min. at 585°C appeared to yield the best stars. Fairly clean 
near field output pattern Is observed when light from a fiber Is Inputted 
Into any of the entrance ports. With light from a fiber (^200 pW) Inputted 
Into one of the entrance ports, the maximum sum total power Intercepted by a 
fiber placed In sequence at the eight exit port was 24 uW (-9.2 dB). There 
are still considerable variations In the outputs of the eight exit ports, 
varying at best by a factor of three. The variations may be due to Imperfection 
scattering, graded index structure boundaries, and/or lack of sufficient mixing 
lengths for the graded Index structure. The cause of these variations is 
presently under investigation. 

We attached an eight-fiber ribbon to each end of a transmission star 
(ion exchange process, 585°C for 20 minutes) with epoxy as described previously. 
During the attachment process, while the epoxy was still liquid, we have 
observed an increase in the throughput to as much as 28 pW (-8.5 dB). Part 
of this Increase can be assigned to the elimination of the air gap; as the 
epoxy hardens, however, a slight decrease (^5 to 10%) of throughput power 
is usually observed. Whether the decrease results from refractive Index 
change of the hardened epoxy or from slight misalignment of the rlbbon-to- 
channel connection as the epoxy hardens Is not known at this time. 
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The 8 by 8 matrices shown In Table la represent the distribution of light 
from each Input fiber to each of the output fibers for the ribbon-star-ribbon 
assembly prior to assembly In a breadboard box. The channels are numbered 
according to Figure 11. Table la shows the initial distribution matrix; 

Table lb shows the distribution where the input and output are interchanged. 
The maximum total throughput of Table la is again 2.4 uW (-9.2 dB), where 
the power in the input fiber is ^200 yW. The average total throughput is 
18.1 yW (-10.4 dB). We note that one of the output channels is much lower 
than the average distribution, particularly with input into channels N and 
#8. We do not know the reason for this behavior at this time. It appears 
that a modification of the outer boundaries of the star, to compensate for 
the graded index of the mixing region, may correct the low output of the 
channel. The computer program for the graded index star is expected to 
clarify the result. 

The lower matrix (with interchanged input and output) shows total outputs 
lower (^-1 to -2 dB) than the upper matrix. At this research stage, we do not 
believe that this is significant. The stars were ion-exchanged in a vertical 
oven with possible thermal gradients in the melt. Inhomogeneous ion exchange 
could have resulted. Future stars are expected to show uniform reciprocal 
distributions as the process is improved. 

The 8x8 distribution matrices of the planar star coupler assembled in a 

i 

box with sixteen fiber couplers and sixteen pigtails as shown in Fig. 7 
was measured and is given in Table Ila, and with Input and output reversed 
in Table lib. We recall that the matrix of Table la is that for the star 
with ribbon pigtails, that of Table Ila, b the same star plus ribbon, but with 
eight input and eight output fiber connectors. The average loss due to the 
two connectors is -'-2.3 dB, or a loss of -1.2 dB per connector. 
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TABLE la. DISTRIBUTION MATRIX FOR 8x8 PLANAR TRANSMISSION STAR #1, 
WITH 2 SETS OF FIBER RIBBONS 


OUTPU' 


INPUTS 
#1 end 


INPUT 
#2 end 
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1 
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2 
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4 
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6 

7 

8 

Total P< 
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2.46 

1.60 

2.26 

1.01 

3.78 

5.99 

5.14 

0.57 

22.31 

3.30 

3.26 

2.48 

5.41 

1.59 

4.05 

2.68 

1.45 

24.22 

0.60 

0.92 

1.65 

3.22 

4.30 

2.43 

2.33 

1.79 

17.24 

0.40 

0.46 

0.58 

0.97 

2.72 

4.56 

3.23 

3.49 

16.41 

2.53 

1.90 

2.02 

1.15 

1.68 

1.80 

1.75 

2.66 

15.49 

1.70 

2.03 

1.98 

1.12 

1.76 

1.67 

2.46 

2.79 

15.51 

6.87 

3.06 

1.38 

1.33 

1.20 

1.23 

0.90 

0.97 

16.94 

1 

0.70 

6.90 

3.05 

1.30 

1.14 

0.93 

0.96 

1.34 

! 16.32 
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18.1 f. 

b. SAME AS la. 

EXCEPT 

INPUT AND OUTPUT ARE REVERSED. 
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1 

I 

2 

3 

4 

5 

6 
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; 3.59 

2.38 

0.85 

0.54 

1.48 

1.24 

1.77 

0.63 

; 12.55 

i 3.29 

i 

3.10 

0.06 

0.69 

1.00 

1.63 

0.85 

1.45 

13.0 

| 3.20 

3.61 

1.42 

0.37 

0.52 

0.78 

0.56 

0.67 

11.13 

! 2.50 

i 

4.60 

3.25 

0.64 

0.30 

0.51 

0.63 

0.69 

1 13.12 

1 

'10.70 

2.50 

1.81 

1.38 

0.55 

0.82 

0.40 

0.46 

! 18.62 

} 

| 3.49 

5.00 

1.43 

1.42 

0.72 

0.88 

0.41 

0.36 

! 13.71 

i 

! 2.37 

2.33 

1.49 

0.56 

1.28 

1.75 

0.47 

0.33 

1 10.58 

j 0.35 

2.30 

0.98 

2.73 

2.51 

1.39 

0.72 

0.50 

j 11.53 








TABLE Ila. DISTRIBUTION MATRIX FOR STAR #1, BREADBOARDED IN A BOX, 
WITH FIBER CONNECTORS 


output 
(white dot 


end) 

1 

2 

3 

4 

5 

6 

7 

8 

Total yW 
Sum Measurement! 

input 

1 

1.6 

1.3 

1.6 

1.3 

2.5 

2.5 

2.8 

.50 

14.1 

\ 

14.40 

2 

1.6 

1.5 

1.5 

3.7 

1.1 

2.9 

1.6 

1.0 

14.9 

15.60 L 

3 

.52 

1.1 

2.0 

3.2 

3.2 

1.5 

1.2 

1.2 

13.9 

13.23 

4 

.33 

.33 

.46 

.58 

1.8 

1.9 

1.5 

1.8 

8.7 

8.70 

5 

1.6 

1.3 

1.0 

.91 

1.3 

1.4 

1.3 

1.8 

10.6 

10.93 i 

6 

1.0 

1.1 

.80 

.56 

.91 

1.6 

.96 

1.6 

8.5 

9.87 

7 

2.0 

1.7 

.62 

.91 

.60 

.41 

.42 

.45 

7.1 

7.65 

8 

.50 

2.7 

1.2 

.60 

.77 

.42 

.50 

.80 

7.5 

8.20 








Avg. 

10.7 
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TABLE lib. SAME AS Ila EXCEPT INPUT AND OUTPUT REVERSED 
output 


Input 
(white 
dot end) 

1 

2 

3 

4 

5 

6 

7 8 

Total mW 
Sum Measurement 

1 

1.6 

1.6 

.57 

.36 

2.0 

.67 

2.0 .73 

9.5 

10.36 

2 

.85 

1.4 

.70 

.26 

.33 

.60 

.54 .67 

5.4 

6.50 

3 

2.1 

1.9 

i.a 

.42 

.59 

.65 

.2? 1.0 

8.2 

8.93 

4 

.72 

.9 

1.7 

.26 

.47 

.40 

.32 .52 

6.3 

6.42 

5 

2.7 

1.5 

2.1 

.81 

.53 

.57 

.24 .40 

8.9 

9.45 

6 

1.2 

2.6 

.90 

1.2 

.86 

.58 

.22 .20 

7.8 

8.50 

7 

.93 

6.1 

6.3 

.90 

1.3 

.68 

.15 .22 

6.6 

8.20 

8 

.13 

1.4 

.93 

2.3 

1.74 

,47 

.38 .30 

7.7 

8.02 








Avg. 

7.6 

8.5 


We wish to note that the matrix measurements of the boxed star were 
sensitive to movement of certain of the Input pigtails and to launch conditions 
of the laser light. Whether this problem arises from distortion at the con- 
nector Is not clear at this time, but the effect will be examined more fully 
and corrected. 

We attached fiber ribbons to another 8x8 transmission star (In exchange 
process, 585°C for 20 minutes) and measured its distribution matrices as 
shown in Tables Ilia and b. The device Is denoted as Star #2. Stars #1 and 
#2 were Ion exchange processed at the same time on the same plate glass sub- 
strate. Comparing Tables I and III shows substantially the same results. One 
difference is the large relative output of channel 8 when channel 1 is excited 
in Table Ilia. Whether this arises from local waveguide defect is not clear 
at this juncture. 

The highest throughput was down *^-96 dB relative to the power in the 
input fiber. We can partition the loss as follows: 

Measured Measured Estimated Estimated 

Insertion = Coupling + Propagation + Internal 
Loss Loss Loss Loss 

or 

9 dB = 4.2 dB + 

The measured coupling loss of 4.2 dB was determined from the experiments 
with parallel channels described above. The remainder of the insertion loss 
is partitioned to propagation loss and the Internal loss of the star struc- 
ture. If the propagation loss is taken to be very small (0.1 dB/cm), an 
estimated Internal loss of 4.0 dB remains. The latter loss may be com- 
pared with the packing fraction loss of 3.8 dB for a linear array of fiber. 

We expect that the fiber to star coupling loss can be reduced to ^0 dB 
with proper matching of index profiles, thereby resulting in an insertion 
loss of ^5 dB. 


(o - 1 ls)( 8 cm ) 


+ [4.0 dB] 
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TABLE Ilia. DISTRIBUTION MATRIX FOR 8x8 PLANAR TRANSMISSION STAR #2 
WITH 2 SETS OF FIBER RIBBONS 

output 


Input 

(X) 

1 

2 

3 

4 

5 

6 

7 

8 

Total yW 
Sum Measurement 

1 

1.1 

.83 

2.5 

2.1 

2.3 

3.2 

1.0 

6.4 

19.4 

21.0 

2 

2.1 

1.7 

3.8 

5.4 

1.7 

3.1 

1.2 

1.7 

20.7 

22.3 

3 

.80 

.95 

2.5 

3.4 

5.2 

2.4 

3.8 

1.2 

20.3 

22.3 

4 

1.4 

1.4 

1.8 

2.4 

2.0 

4.7 

2.5 

3.4 

19.6 

17.8 

5 

.87 

1.4 

1.3 

1.2 

1.1 

2.4 

1.4 

3.3 

13.0 

12.3 

6 

6.1 

2.6 

1.6 

1.6 

2.0 

1.3 

1.4 

.98 

17.6 

16.3 

7 

7.4 

5.6 

2.1 

.76 

1.2 

.70 

.63 

.76 19.2 

18.8 

8 

.68 

8.0 

3.3 

3.6 

1.7 

.96 

1.4 

1.7 

21.3 

21.2 








Avg. 


18.9 

19.0 


TABLE Illb. SAME AS Ilia EXCEPT INPUT AND OUTPUT ARE REVERSED 
output 

(X) Total yW 


i nput 

1 

2 

3 

4 

5 

6 

7 

8 

Sum 

Measurement 

1 

r>» 

oo 

♦ 

1.1 

.42 

.75 

.61 

3.8 

13.4 

.78 

21.7 

22.4 

2 

.48 

.82 

.46 

1.3 

.77 

1.5 

5.6 

7.0 

17.9 

17.6 

3 

1.7 

1.9 

2.3 

1.8 

.93 

1.4 

2.7 

2.3 

15.0 

14.6 

4 

1.7 

1.1 

1.6 

2.4 

1.2 

1.4 

.84 

2.7 

12.9 

13.5 

5 

2.1 

1.7 

3.3 

3.5 

.95 

1.5 

.82 

1.8 

15,7 

15.7 

6 

2.1 

6.8 

2.2 

4.0 

.96 

.74 

.62 

.99 

18.4 

18.6 

7 

3.8 

3.0 

3.3 

1.5 

.78 

.99 

.42 

.77 

14.6 

15.5 

8 

4.6 

2.3 

1.4 

2.8 

1.9 

.75 

.61 

1.1 

15.5 

14.6 








Avg. 


16.4 

16.6 
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CONCLUSIONS 

• Manufacture of planar distribution couplers for optical fibers 
where the planar structure Is fabricated by an ion diffusion 
process has been demonstrated. 

• Mechanically solid coupler-fiber ribbon-connector assembly 
possible for future spacecraft applications have been demonstrated. 

• No restriction on the number of fiber terminals appears to be 
valid. Multiple couplers may be used in series. 

t Further improvement of the couplers to -5 dB transmission loss 
[plus connector loss (M.2 dB/connector)] appears to be attainable. 

• Precise equal partitioning of light appears to be possible. 
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DESIGN OF PLANAR STAR COUPLERS 
FOR FIBER OPTIC SYSTEMS 

0. G. Rawer 

RECEDING PAGE BLANK HOT FU*» 

ABSTRACT 

A horn structure within a planar fiber optic star 
coupler has been shown by computer design to 
significantly reduce If not eliminate the packing 
fraction loss associated with the coupler geometry. 

The ray tracing design procedure is used to 
characterize the partitioning of light to the various 
ports and the losses are tabulated. An eight port 
planar star is used to illustrate the design 
capabilities. 

INTR ODUCTION 

Low loss multimode fiber optic components are a key to successful 
multimode fiber optic systems. The technique of fusing individual fibers 

1 O 1 

together '* to make low loss T couplers has resulted in reproducible 
couplers but the yield of packaged devices is low, suggesting high unit 
costs. Fiber optic star couplers^ are successfully made by a number of 
techniques, however the packing fraction loss remains a problem. A 
planar approach of making these couplers is being pursued. It Is believed 
that high yield and low cost can be achieved. This paper describes a 
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0 

ray-tracing Investigation of the multimode planar star coupler Illustrated jj 

L; 

In Fig. 1. This Investigation was made In order to characterize the 

partitioning of light Into the planar branches as a function of Index |j 

profile, device dimensions and configuration; the losses are also 

) 

tabulated. Computer ray tracing has been used to Identify the Important 

design parameters for the planar star. Example parameters are the ) \ 

required mixing length for equal partitioning of light Into each of the 

i .* 

channels and the horn angle required to significantly reduce the packing [j 

fraction losses. It should be pointed out that straight line ray tracing 
requires that the planar guides and the substrate be uniform in refractive 
Indices. The ray tracing analysis is limited to such planar structures. 

RAY TRACING PROCEDURE 

I ! 

( j 

In order to accurately simulate the devices, the rays must statistically 
simulate the total of the modes In the guides. Since the structure that Is \ 

analyzed is a step Index configuration, we have elected to choose the rays 
randomly such that the angular extent of the numerical aperture Is filled 
uniformly. Although It is also possible to randomly choose the Initial x 
and y position of the rays In some z plane, we elected to divide the input 
plane into a rectangular mesh and trace a single ray from the center of each j 

rectangle. Using this algorithm and starting the random number generator 
at different positions does not affect the final result significantly 
(<1% variation) If a large number (>5000) of rays are traced. That Is, 

‘ f 

some average properties associated with the motion of a large number of 
Identical systems differing over a range of Initial conditions Is obtained. 
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Given the initial values for position (x,y) and direction cosines 
(k, 1, m) of a ray at some 2 Q {*0) it Is possible, in general, to solve 
for its position and direction cosines at any other value of z. Although 
ray tracing Is well understood the procedure Is described. 

Snell's law describes the refraction of light through an Interface 


(Fig. 

1) and is given by 

a 


or 

n^ Sg x E ■ nj x G 

(1) 

where 

(n^ “ x (a » 0 

(2) 


Is a unit vector in thd direction of the ray impinging on 
the interface 

G is a unit vector in the direction of the local normal to 
the interface 

n^ qnd n^ are the refractive indices on the two sides of the 
Interface 

$2 is a unit vector in the direction of the ray after refraction 
and/or reflection (n-j * n 2 ) at the interface. 

Equation (2) is a simple vector relation which can be solved exactly; it is 

satisfied if the vector in the parentheses is col inear with G, thus 

allowing one to write 


n£^2 ” n^i = TG 

(3) 

n 2^2 * n l^l * 

(4) 


where r is yet an unknown quantity. The dot product of Eq. (4) with Itself 
eliminates the unknown and gives the following quadratic equations in T; 
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or 


r 2 - 2r(n^ • IS) + (rtj 2 - n 2 2 ) B 0 


(5) 


r * n^Sj • G 1 t /(n 1 f 1 S) z - n^ + n 2 2 . (6) 

By the appropriate choice of r, all situations that are encountered at a 
refractive Interface can be accurately described. If the square root is 
Imaginary (n-j t n 2 ), the ray Is totally Internally reflected and the 
positive sign Is chosen after setting n 1 * n 2 

^2 * ^1 " ^ 2S 1 * (7) 

If the square root Is real, the smaller absolute value of r will give the 
S 2 corresponding to the physical situation of refraction (Fig. 3). 
Substitution of this value for r Into Eq. (4) gives the new ray direction 
? 2 . In the analysis rays of this sort are lost from the device and thus 
summed to evaluate the losses. 

The algorithm for ray tracing through a device Is 

a) initialize ray, 

b) determine surface of intersection, 

c) apply Eq. (7) If square root In Eq. (6) Is Imaginary. 

If square root is real, sum the ray as lost and go 
to (a) 

d) repeat (b) and (c) until ray exits device, 

e) sum ray as exiting appropriate port and to to (a). 

The parameters of Interest for star design are (see Fig. 1): 

$ » half angle of taper angle, 

Xy - half width of wall, 

Xg * half width of channel, 

Yg * half thickness of film, 
z * mixing length, 

• 

Numerical aperture of the guide, 

Index of the surrounding region. 

Number of ports N. 
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In the above algorithm the ray intersections with the guide boundaries 
remains to be defined. There are five Independent sections (input guides, 
output guides. Input horns, output horns, and mixing region) of the planar 
star. Each has to be handled separately with the appropriate coordinate 
transformation between regions. The coordinate system for each Individual 
guide (horn, mixing region, etc.) Is a 2 axis that bisects the guide with 
transverse x, and y axis. The x axis Is parallel to the plane of the substrate. 

Let the ray position be given by 

F^xT + yJ+zF (8) 


and the ray direction by 

T e lT + mT+nF (9) 

\ 

If we are in any of the three^ uniform regions (Input guide, output guide 
or the mixing region) the next boundary intersection, Pj, is given by 

P j * P + T 0 (10) 

where D is the smallest of the following three distances: 
a) the distance to the x plane 


OX « (x -tXg/|t|)/t 
b) the distance to the y plane 

DY * (y - m Yq/ |**»J )/m 

and c) the distance to the end of the guide 

OZ • (z - z rax )/n 

where 2x 6 , 2y G and z {nax are the guide dimensions. For example. In the 
mixing region 

x 6 - (X B + x t )n 
y g ■ y b 

z max ’ z mfx 


(ID 

( 12 ) 

(13) 


(14) 

(15) 

(16) 
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In the horns the next boundary Intersection (Eq. 10) Is given 
by the shortest distance between the y plane (Eq. 12), the end of the 
horn (Eq. 13) z^ ■ Xy/tan 0 or one of the two horn boundaries. 


D - 


X - (Mz + B ) 
Mn * £ 


where for the Input horns 

M * t tan 0 

B * - Xg 

and for the output horns 

M « - tan(-e) 

B * t (Xg + Xj.) 


This completes the ray tracing algorithm. 


07 ) 

( 18 ) 

09 ) 

(20) 

(21) 


RAY TRACE ANALYSIS 

n the ray trace program, rays can be Input at any of the N ports; 
the program outputs are the number of rays that exit each of the ports 
and the number of rays that are lost. Figure 3 shows the result of varying 
the mixing length, Z, for an eight-port star coupler; each line represents 
the number of rays that exit a given port. The number of rays to each of 
the ports with ports 1 through 4 for Input and Z/Xg * 2750 are given in 
Table 1. The number of rays lost or not satisfying the condition for total 
Internal reflection at ar\y one of the guide wall Intersections Is also 
given In Table 1. A total of 10,000 rays were traced for each Input. 

The major loss In conventional star couplers Is the packing fraction 
loss. The purpose of the horn structure In the planar star coupler Is to 
provide a means to reduce this loss. The losses for an eight port star 
have been calculated with the horn angle, 6, as a variable. The result 
of this calculation Is shown In Fig. 4. For a horn angle equal to or 
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TABLE 1. 


“Input 

port 

^^Dutput 

port 

1 

2 

3 

4 

1 

1195 

999 

927 

905 

2 

967 

942 

922 

957 

3 

948 

934 

964 

1065 

4 

952 

1014 

1072 

1053 

5 

891 

967 

1084 

1077 

6 

966 

991 

992 

1082 

7 

1020 

988 

1002 

918 

8 

1162 

1107 

997 

884 

Loss 

1999 

2058 

2040 

2059 

Total 

10,000 

10,000 

10,000 

10,000 


greater than half the numerical aperture (NA/2) of the guide the loss Is 
equal to the structure packing fraction. In this case the less Is 

tyUy + X B ) » .75 (21) 

At reduced angles the loss Is substantially reduced. The packing fraction 
loss for horn angles larger than NA/2 can be easily understood. Any ray 
that strikes the horn boundaries will have Its direction changed by an 
angle equal to twice the horn angle. Hence for any horn angle larger than 
NA/2 the ray angle after striking the horn will be larger than the guide NA 
and lost. Reduced horn angles result in an adiabatic transition and 
reduced loss. 
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In conclusion the planar star coupler design presented shows that 
the packing fraction loss of the planar star copuler can be decreased. 

This is achieved by tapering the Input guides with a taper angle less 
than half the numerical aperture of the guide. The mixing region length 
can be chosen for equal coupling to all ports. An experimental Investigation 
is presently underway and the results will be reported In a separate article. 

The author appreciates the critical reading of the manuscript and 
suggestions of Dr. C. K. Asawa. The stimulating discussions of Dr. M. K. 
Barnoskl and Dr. G. L. Tangonan are also acknowledged. 
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FIGURE CAPTIONS 


1. Schematic of a transmissive planar star coupler. Rays from 
the left are partitioned in the mixing region and exit from 
the guides at the right. The F. indicate the fractional 
transmission to the i output port. 

2 . Illustration showing that the appropriate choice of r in the 
solution of Snell’s law is 

3. Fractional transmission of an eight-port planar star coupler 
versus the mixing length. The fractional transmission to port i 
is labeled F^. 

4. Fractional loss of an eight port star versus the horn angle, 0 . 
The loss of 75% is the structure packing fraction. 





FRACTIONAL LOSS 
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Ion exchange in a melt of lithium and potassium sulfates is utilized to form couplers for multimode fiber optics Flic 
results of coupling measurements on asymmetric Y -branch couplers are given. 


Foi multimode fiber optic systems couplers or taps 
are necessary for data distribution to several terminal 
users. Couplers are presently made with the help of 
optical fibers (fusing, tapering, lapping and gluing) or 
with microoptic components such as microlenses and 
beam splitters. Good couplers have been demonstrated 
with all these methods. 

Multimode couplers can also be made by planar 
processing. In principle, the planar approach utilizing 
photolithographic processing is promising because of 
its versatility and high reproducibility. Auracher et al. 

1 1 1 have reported the first planar branching networks 
formed in photo-polymer materials of 100 pm thick- 
ness. We have formed planar coupling structures for 
multimode guides using the ion exchange process de- 
scribed by Chartier et al. [2]. In this process an eutec- 
tic mixture of Li 2 S0 4 and K0SO4 is heated under an 
oxygen atmosphere to 580°C. A sodium glass slide is 
suspended over the melt for 30 minutes to reach ther- 
mal equilibriation with the melt. It is next dipped in- 
to the melt for 20 minutes, then again suspended 
over the melt for 10 minutes to avoid any thermal 
shock. Planar waveguides. 100 pm deep, are made by 
this process. The coupler structures are formed by 
masking with a thick Al film ( 1 to 2 pm). The Al 
mask is subsequently removed by dipping the slide in 
hot 6 M HCI solution. Degradation due to aging was 
not observed, unlike the Ag ion exchange processes 
in glass [3). 

An investigation of reactions between the melt and 
the glass during the process was carried out. Our study 
on the melt showed that SO3 was continually produced. 
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U2SO4, which is less thermally stable than k-.S() 4 . 
was decomposing to Li 2 0. which is corrosive to S1O1. 
The addition of Li 2 0 to Si0 2 to form L.i 2 SiO; ex- 
plains the observed thinning of the slide by 2 pm alter 
20 minutes of exchange. It would be difficult to at- 
tribute this thinning solely to a simple exchange pto- 
cess. 

Whether the increase in refractive index is due ex- 
clusively to an ion exchange of Li 4 for Na + has not 



300, m SPAC NO 
Bt TWEEN 
GUIDES 
1 7!> . 111 WIDE 



OUTPUTS OF 1 COUPLER 
1300 uni CENTER TO CENTER) 


Fig. I . Planar coupler array formed by inn exchange in a glaw 
slide. 
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been ascertained. Other exchange processes such as 
K + for Na +(4 * or addition of L^Q [5 J may be occurr* 
ing. The relative contributions of these other mecha- 
nisms have not been determined. Experiments using 
selected lithium and potassium salts (LiCl and KC1, 
for instance) are being conducted to determine the 
importance of these effects. 

In fig. 1 we show a top view of an ion exchange 
slide with coupling structures formed using 1 , 2, 3, 4 
and 5° branches. The coupling branches are continued 
parallel to the main branch when 300 pm separation is 
attained. This facilitates the power measurement be- 
cause the optical fibers placed at the output ports may 
be perpendicular to the slide end face. The output of 
the 1° branch is shown also in fig. 1. Using graded in- 
dex fibers (Corning) of 65 pm core size the coupling 
ratios were measured from fiber input to fiber output 
with the planar coupler in between. Mode strippers 
were used to ensure that only guided modes were 
measured. The results are shown in fig. 2, where for 
two outputs P x and P 0 we have graphed P\!(P\ + Pq) 
where P j is the branch output and P 0 is the straight 
output. A clear trend of decreasing coupling with in- 
creased branching angle is observed. The error bars 
are due mainly to the experimental difficulties in 
fiber placement precision and processing variations. 

To demonstrate further the versatility of the planar 
coupler approach a tapered coupler structure has been 
formed as shown in fig. 3. The coupling is controlled 
by varying the parameters ® 2 a °d ®3* P° r ^i 



Fig. 2. Coupling ratio results for different branch angles. 

* 75 , <#2 * 50 ,d 3 * 100pm and 6 X = 2 °,0 2 = 0.66°, 
and = l°,a coupler with 20: 1 tap ratio was mea- 
sured. 

The ion exchange process seems well suited for 
fiber optic coupler applications, in fact, the guide los- 
ses we have observed are quite low. We estimate that 
the guide losses are 0.1 to 0.5 dB/cm depending on 



Fig. 3. Tapered couplers of planar variety formed by ion exchange. Coupling is due predominantly to radiation mode excitation 
by the taper. 
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processing variations. The dominant scattering mecha- 
nism is due to surface cracks, which have been largely 
eliminated by preprocessing of the slide. The crack 
density is diminished by a process involving a 5 min 
etching of the slide in 5% aqueous HF. In addition, a 
fire polishing of the slides prior to ion exchange re- 
duces the final crack density. A detailed description 
of these results will be published later [6] . 

Several problem areas must be addressed which 
deal with optimization of the coupling efficiency from 
fiber to planar coupler, device packaging (affixing the 
pigtails), and guide formation reproducibility. Our 
best results on excitation of the channel guides with 
graded index fibers give only a 65% coupling efficiency. 
The processing used was 30 min diffusion at 580°C 
with a 65 um channel width opening. At present this 
is the dominant loss figure. It is clear from fig. I that 
the output intensity profile must be modified to bettei 
approximate that of a fiber to improve the coupling. 
Careful control of the diffusion conditions will cer- 
tainly enhance this figure, especially when buried 
structures are obtained by a second exchange process 
wherein K is re-introduced into the glass, subsequent 
to the Li exchange. 


In conclusion we have demonstrated the feasibility 
of the ion exchange process in forming deep multi- 
mode guide structures for fiber optic coupler applica- 
tions. 

The authors wish to thank Drs. M.K. Bar noski. 

O.G. Ramer, C.K. Asawa, C.R. Blair and R.C. Pastor 
for their encouragement and advice. 
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